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ABSTRACT 

A  X  ^  analysis  of  standard  model  accretion  disk  synthetic  spectrum  fits  to  combined  Far  Ultraviolet  Spectroscopic 
Explorer  and  Space  Telescope  Imaging  Spectrograph  spectra  of  V3885  Sagittarius,  on  an  absolute  fiux  basis,  selects 
a  model  that  accurately  represents  the  observed  spectral  energy  distribution.  Calculation  of  the  synthetic  spectrum 
requires  the  following  system  parameters.  The  cataclysmic  variable  secondary  star  period-mass  relation  calibrated 
by  Knigge  in  2006  and  2007  sets  the  secondary  component  mass.  A  mean  white  dwarf  (WD)  mass  from  the  same 
study,  which  is  consistent  with  an  observationally  determined  mass  ratio,  sets  the  adopted  WD  mass  of  0.7  Mq,  and 
the  WD  radius  follows  from  standard  theoretical  models.  The  adopted  inclination,  i  =  65°,  is  a  literature  consensus, 
and  is  subsequently  supported  by  x^  analysis.  The  mass  transfer  rate  is  the  remaining  parameter  to  set  the  accretion 
disk  Teff  profile,  and  the  Hipparcos  parallax  constrains  that  parameter  to  M  =  (5.0  zb  2.0)  x  10“^  Mq  yr“^  by  a 
comparison  with  observed  spectra.  The  fit  to  the  observed  spectra  adopts  the  contribution  of  a  57,000  ±  5000  K 
WD.  The  model  thus  provides  realistic  constraints  on  M  and  Tgff  for  a  large  M  system  above  the  period  gap. 

Key  words:  novae,  cataclysmic  variables  -  stars:  individual  (V3885  Sagittarius)  -  white  dwarfs 
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1.  INTRODUCTION 

Cataclysmic  variables  (CVs)  are  semi-detached  binary  stars 
in  which  a  late  main-sequence  star  loses  mass  onto  a  white  dwarf 
(WD)  by  Roche  lobe  overflow  (Warner  1995).  In  non-magnetic 
systems,  the  mass  transfer  stream  produces  an  accretion  disk 
with  mass  transport  inward  and  angular  momentum  transport 
outward,  driven  by  viscous  processes.  The  accretion  disk  may 
extend  inward  to  the  WD;  the  outer  boundary  extends  to 
a  tidal  cutoff  limit  imposed  by  the  secondary  star  in  the 
steady  state  case.  If  the  mass  transfer  rate  is  below  a  certain 
limit,  the  accretion  disk  is  unstable  and  undergoes  brightness 
cycles  (outbursts),  and  if  above  the  limit,  the  accretion  disk 
is  stable  against  outbursts.  The  latter  objects  (which  have  no 
recorded  outburst  of  any  type)  are  called  nova-like  (NL)  systems. 
As  shown  by  current  studies  (Ribeiro  &  Diaz  2007;  Hartley 
et  al.  2005,  2002;  Woods  et  al.  1992;  Metz  1989;  Haug  & 
Dreschel  1985;  Guinan  &  Sion  1982;  Cowley  &  Crampton 
1977),  V3885  Sgr  is  a  NL  system  (see  Warner  1995  for  a 
detailed  discussion).  V3885  Sgr  has  a  Hipparcos  (Perryman 
et  al.  1997)  parallax  of  9.11  ib  1.95  mas.  The  orbital  period 
is  0^20716071  (Ribeiro  &  Diaz  2007),  and  the  inclination  is 
i  ~  65°  (Haug  &  Dreschel  1985;  Hartley  et  al.  2005;  Ribeiro  & 
Diaz  2007). 

NL  systems  are  of  special  interest  because  they  are  expected 
to  have  an  accretion  disk  radial  temperature  profile  given  by 
an  analytic  expression  (Erank  et  al.  1992,  Equation  (5.41); 
hereafter  EKR)  which  defines  the  so-called  standard  model; 
the  expression  includes  the  mass  transfer  rate  M  as  an  explicit 
variable.  In  NL  systems,  which  are  above  the  period  gap 

*  Based  on  observations  made  with  the  NASA/ESA  Hubble  Space  Telescope, 
obtained  at  the  Space  Telescope  Science  Institute,  which  is  operated  by  the 
Association  of  Universities  for  Research  in  Astronomy,  Inc.  under  NASA 
contract  NAS5-26555,  and  the  NASA-CNES-CSA  Far  Ultraviolet  Explorer, 
which  is  operated  for  NASA  by  the  Johns  Hopkins  University  under  NASA 
contract  NAS5-32985. 


(Howell  et  al.  2001),  the  accretion  disk  dominates  the  system 
spectrum.  Eitting  a  synthetic  spectrum,  based  on  the  analytic 
model,  a  proxy  of  the  accretion  disk  temperature  profile  to 
an  observed  spectrum  potentially  determines  M,  a  physical 
parameter  of  basic  importance,  since  it  controls  the  evolution  of 
CV  systems  (Howell  et  al.  2001).  But  the  analytic  expression 
also  is  an  explicit  function  of  the  WD  mass,  M^d,  and  it  must 
be  determined  independently. 

By  the  basic  paradigm  of  CV  evolution  (Howell  et  al.  2001), 
there  is  a  close  relation  between  orbital  period  and  donor 
mass.  In  a  recent  paper,  Knigge  (2006,  2007)  has  calibrated 
the  relations  among  P,  M2,  Teff,25  and  donor  spectral  type. 
Thus,  given  an  orbital  period  below  the  tabular  upper  limit  of 
6  hr,  M2  can  be  determined  independently  of  other  parameters. 
If  there  is  an  observationally  determined  mass  ratio,  the  WD 
mass  follows  directly. 

A  normalizing  factor  is  necessary  in  comparing  synthetic 
spectra  with  observed  spectra,  and  the  factor  depends  on  the 
distance  to  the  system.  If  the  distance  is  known  independently, 
an  absolute  determination  of  M  potentially  is  possible. 

A  complication  in  fitting  a  synthetic  spectrum  to  an  observed 
spectrum  is  due  to  the  fact  that  the  WD  may  contribute  detectably 
to  the  system  spectrum;  evaluating  the  WD  contribution  requires 
the  WD  Teff.  A  limited  number  of  CV  systems  have  observa¬ 
tionally  determined  values  for  their  WD  Tgff  (Sion  1999;  Winter 
&  Sion  2003),  with  extremely  few  being  NLs.  Because  of  the 
large  M  in  NL  systems,  the  accretion  disk  dominates  the  system 
spectrum  in  the  UV  and  optical,  and  makes  determination  of  the 
WD  Teff  difficult.  Thus,  the  Teff  of  an  additional  NL  WD  is  of 
obvious  interest. 

2.  THE  FUSE  AND  STIS  SPECTRA 

Table  1  presents  the  observations  log.  The  Space  Telescope 
Imaging  Spectrograph  (STIS)  spectra  are  the  same  as  three  STIS 
spectra  described  by  Hartley  et  al.  (2002). 
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Table  1 

FUV  Observations  of  V3885  Sgr:  HST  &  FUSE  Spectra 


Date 

Telescope/ 

Exp  Time 

Data  Set 

Filter/Grating/ 

Operation 

Calibration 

(yyyy/mm/dd) 

Instrument 

(s) 

Aperture 

Mode 

Software 

2000  May  24 

FUSE 

12392 

P1870101000 

LWRS 

TTAG 

CalEUSE  3.0.7 

2000  Apr  30 

STIS 

2195 

05BI04010 

E140M/0.2  X  0.2 

TIME-TAG 

CALSTIS  2.22 

2.7.  The  FUSE  Spectrum 

The  Far  Ultraviolet  Spectroscopic  Explorer  {FUSE)  is  a  low- 
earth  orbit  satellite,  launched  in  1999  June.  Its  optical  system 
consists  of  four  optical  telescopes  (mirrors),  each  separately 
connected  to  a  different  Rowland  spectrograph.  The  four  diffrac¬ 
tion  gratings  of  the  four  Rowland  spectrographs  produce  four 
independent  spectra  on  two  microchannel  plates  (detectors). 
Two  mirrors  and  two  gratings  are  coated  with  SiC  to  provide 
wavelength  coverage  below  1020  A,  while  the  other  two  mir¬ 
rors  and  gratings  are  coated  with  A1  and  LiE.  The  Al-fLiE  coat¬ 
ing  provides  about  twice  the  reflectivity  of  SiC  at  wavelengths 
>1050  A,  and  very  little  reflectivity  below  1020  A.  These  are 
known  as  the  SiCl,  SiC2,  LiEl,  and  LiE2  channels  producing 
a  total  of  8  (partially  overlapping)  spectral  segments:  SiC  la, 
SiC  lb,  SiC2a,  SiC2b,  LiEl  a,  LiElb,  LiE2a,  and  LiE2b  for  each 
exposure. 

Eour  FUSE  observations  of  V3885  Sgr  prospectively  are 
available  for  analysis.  An  observation  on  2002  August  10  has  a 
short  exposure  time  and  has  an  inferior  S /N  ratio  compared  to 
other  exposures.  An  exposure  on  2000  May  24  (UT:01:51:55) 
has  the  longest  exposure  time  and  occurs  within  25  days  of  our 
chosen  STIS  spectrum.  We  chose  this  spectrum  for  analysis. 
The  exposure  used  the  30''  x  30"  LWRS  Large  Square  Aperture 
in  TIME  TAG  mode.  The  telescope  collected  11  exposures 
obtained  during  1 1  consecutive  orbits  of  the  FUSE  spacecraft, 
totaling  12865  s  of  raw  exposure  time  for  the  88  spectral 
segments;  we  downloaded  the  spectra  from  the  MAST  archive 
and  processed  them  as  described  below. 

The  data  were  processed  with  the  latest  and  flnal  version  of 
CalEUSE  (v3.2.1;  Dixon  et  al.  2007).  The  main  change  from 
previous  versions  of  CalEUSE  is  that  the  data  are  maintained 
as  a  photon  list  (the  intermediate  data  file:  IDE)  throughout  the 
pipeline.  Bad  photons  are  flagged  but  not  discarded,  so  the  user 
can  examine.  Alter,  and  combine  data  without  re-running  the 
pipeline.  Event  bursts  (short  periods  during  an  exposure  when 
high  count  rates  are  registered  on  one  or  more  detectors)  are 
processed  automatically  in  this  latest  version  of  CalEUSE.  The 
bursts  exhibit  a  complex  pattern  on  the  detector;  their  cause, 
however,  still  is  unknown  (it  has  been  confirmed  (Dixon  et  al. 
2007)  that  they  are  not  detector  effects). 

The  11  exposures  all  have  essentially  identical  continuum 
flux  levels,  and  there  is  no  detectable  change  in  flux  from 
exposure  to  exposure.  The  main  features  that  vary  from  exposure 
to  exposure  are  the  depth,  width  and  shape  of  the  absorption 
lines.  The  next  step  is  to  co-add  the  11  exposures,  for  the  eight 
individual  segments.  The  count  rate  indicates  screening  of  22  s 
of  high  voltage  (HV)  for  the  eight  spectral  segments.  In  addition, 
there  was  a  loss  of  good  exposure  time  due  to  event  bursts  of 
472  s,  409  s,  390  s,  and  511  s  in  segments  la’s,  lb’s,  2a’s,  and 
2b’s,  respectively.  We  obtain  a  total  of  12310  s,  12373  s,  12392  s, 
and  12271  s  for  segments  la’s,  lb’s,  2a’s,  and  2b’s,  respectively, 
including  more  than  8200  s  of  NIGHT  time  (i.e.,  when  the 
spacecraft  is  located  in  the  shadow  of  the  earth).  Since  the 
total  good  collection  time  varies,  we  weighted  the  segments 
accordingly  (see  below). 


During  the  observations,  fine  error  sensor  A  (which  images 
the  LiEl  aperture)  was  used  to  guide  the  telescope.  The  spectral 
regions  covered  by  the  spectral  channels  overlap,  and  these 
overlap  regions  are  then  used  to  renormalize  the  spectra  in  the 
SiCl  (a  &  b),  LiE2  (a  &  b),  and  SiC2  (a  &  b)  channels  to  the  flux 
in  the  LiEl  (a  &  b)  channel.  The  low  sensitivity  portions  (usually 
the  edges)  of  each  channel  are  discarded.  In  addition,  there  exists 
a  narrow  dark  stripe  of  decreased  flux  in  the  spectra  running  in 
the  dispersion  direction,  affectionally  identified  as  the  “worm,” 
which  can  attenuate  as  much  as  50%  of  the  incident  light  in  the 
affected  portions  of  the  spectrum  due  to  shadows  of  the  wires 
in  a  grid  above  the  detector.  Because  of  the  temporal  changes  in 
the  strength  and  position  of  the  worm,  CalEUSE  (Dixon  et  al. 
2007)  cannot  correct  target  fluxes  for  its  presence.  Therefore, 
we  carried  out  a  visual  inspection  of  the  FUSE  channels  to 
locate  the  worm,  and  we  manually  discarded  the  portion  of  the 
spectrum  affected  by  the  worm:  ~1 130-1 189  A  in  the  LiElb 
channel. 

We  combined  the  individual  exposures  and  channels  to 
create  a  time-averaged  spectrum,  weighting  the  flux  in  each 
output  datum  by  the  exposure  time  (as  mentioned  earlier)  and 
sensitivity  of  the  input  exposure  and  the  channel  of  origin.  The 
flnal  product  is  a  spectrum  that  covers  almost  the  full  FUSE 
wavelength  range  905  A-1187  A;  because  we  disregarded  that 
segment  of  the  spectrum  affected  by  the  worm,  the  spectrum 
ends  at  ~  1 1 82  A.  The  absorption  lines  in  the  flnal  spectrum  are 
blueshifted  by  about  0.4  A.  Based  on  the  Ribeiro  &  Diaz  (2007) 
ephemeris,  the  FUSE  spectrum  extends  from  orbital  phase  0.29 
to  orbital  phase  1.01. 

2.2.  The  STIS  spectra 

The  STIS  spectra  (from  the  MAST  archive)  were  taken  in 
TIME-TAG  mode,  with  the  EUV-MAMA  detector.  The  E140M 
echelle  grating  was  used  (with  a  central  wavelength  of  1425  A 
and  the  0.2  x  0.2  aperture).  The  STIS  spectra  were  processed 
with  CALSTIS  version  2.22.  The  calibrations  were  performed 
using  the  STSDAS  package.  Each  STIS  spectrum  includes 
44  echelle  spectra,  each  one  covering  roughly  18  A.  The 
echelle  spectra  overlap  except  in  the  longer  wavelengths.  As 
a  consequence,  when  the  echelle  spectra  are  merged  to  produce 
the  flnal  spectrum,  there  are  four  gaps  in  the  longer  wavelengths. 
The  gaps  are  about  "^0.5  A  wide  and  are  at  intervals  of  about 
1 8  A.  We  have  removed  the  gaps  by  hand. 

The  first  and  third  STIS  spectra  have  closely  similar  continua 
and  differ  slightly  in  the  shape  of  their  absorption  lines,  with  a 
blueshift  of  1.3  A  and  0.6  A.  The  second  STIS  spectrum  has 
a  continuum  flux  level  14%  larger,  and  its  absorption  lines  are 
blueshifted  by  only  0.4  A.  Based  on  the  ephemeris  of  Ribeiro 
&  Diaz  (2007),  the  Table  1  STIS  spectra  extend  in  orbital  phase 
from  0.67  to  0.81,  0.33  to  0.47,  and  0.35  to  0.49,  respectively. 
The  STIS  spectra  cover  the  interval  from  1 150  A  to  1720  A. 

2.3.  The  Combined  Spectra 

The  FUSE  and  first  STIS  spectra  of  V3885  Sgr  are  shown 
in  Eigure  1  before  correction  for  reddening;  the  spectra  are 
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Figure  1.  FUSE  plus  STIS  combined  spectra  with  line  identifications. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


characterized  by  broad  and  deep  absorption  lines  from  the 
source  (disk  plus  wind/chromosphere).  These  broad  absorption 
lines  are  blueshifted  by  about  1.5  A  and  2.5  A  in  the  FUSE 
and  STIS  spectra  respectively;  they  are  about  4  A  in  width 
(~1000  km  s“^)  and  exhibit  the  same  asymmetry  in  both  the 
FUSE  and  STIS  spectra.  The  red  edge  of  the  absorption  is  steep, 
while  the  blue  edge  has  a  more  gradual  slope,  resembling  the 
blueshifted  absorption  part  of  a  P  Cygni  profile.  However,  only 
the  Civ  (1550)  and  Nv  (1240)  lines  exhibit  a  true  P  Cygni 
profile,  while  all  the  other  broad  absorption  lines  do  not.  All 
these  lines  belong  to  high  ionization  species:  Niii,  Niv,  Nv, 
Si  III,  Si  IV,  C  III,  C IV,  O  VI,  P  V,  and  He  ii. 

In  addition  to  the  broad  lines,  there  are  many  narrow  and 
deep  absorption  lines.  They  have  a  width  varying  between  about 
0.15  A  (e.g.  Oi)  and  0.3  A  (for  Hi)  and  are  also  blueshifted, 
but  by  a  much  smaller  amount  (^^0.2  A)  than  the  broad  lines. 
Many  of  the  narrow  lines  belong  to  low-order  ionization  species; 
however,  some  also  belong  to  higher  ionization  species.  Some 
of  the  sharp  lines  are  known  to  be  interstellar,  such  as  Ar  i  and 
A1 II.  For  simplicity,  we  put  all  these  lines  in  the  ISM  category. 
Additional  details  concerning  the  STIS  spectra  are  in  Hartley 


et  al.  (2002).  Tables  2  and  3  list  all  the  lines,  separated  into 
FUSE  and  STIS  lines. 

From  the  AAVSO  Web  site,^  we  determine  that  during  the 
month  of  2000  April  and  through  2000  May,  the  V3885  Sgr 
visible  magnitude  varied  by  only  0.2.  The  first  of  the  Table  1 
STIS  spectra  falls  within  this  time  frame,  and  it  is  closest  in 
time  to  the  FUSE  spectrum.  Also,  the  orbital  phases  of  the  two 
spectra  are  nearly  the  same.  Finally,  the  first  STIS  spectrum  most 
closely  matches  the  FUSE  spectrum  in  the  overlap  region.  All 
of  the  spectra  have  been  corrected  for  interstellar  reddening  of 
E{B  —  V)  =  0.02  (Section  3).  We  choose  the  first  STIS  spectrum 
to  combine  with  the  FUSE  spectrum  for  analysis.  Figure  2  shows 
the  combination  of  the  FUSE  spectrum  with  the  three  STIS 
spectra.  The  upper  panel  shows  the  overlap  region.  The  FUSE 
spectrum  is  red,  the  first  STIS  spectrum  is  cyan,  the  second  STIS 
spectrum  is  magenta,  and  the  third  STIS  spectrum  is  brown.  For 
optimum  fit  to  the  FUSE  spectrum,  the  STIS  spectra  have  been 
divided  by  0.98,  1.12,  and  0.94,  respectively.  The  top  panel 
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Table  2 

FUSE  Lines 


Line 

Identification 

Wavelength 

(A) 

^0 

(A) 

Origin 

Hi 

914.13 

914.29 

Ism 

Hi 

914.43 

914.58 

Ism 

Hi 

914.79 

914.92 

Ism 

Hi 

915.19 

915.33 

Ism 

Hi 

915.68 

915.82 

Ism 

Hi 

916.30 

916.43 

Ism 

Oi 

916.73 

916.82 

Ism 

Hi 

917.05 

917.18 

Ism 

Hi 

918.00 

918.13 

Ism 

Hi 

919.22 

919.35 

Ism 

Hi 

920.84 

920.96 

Ism 

Hi 

921.75 

921.86 

Ism 

Hi 

923.01 

923.15 

Ism 

Niv 

-922 

923.2 

S 

Oi 

924.83 

924.95 

Ism 

Oi 

925.33 

925.45 

Ism 

Hi 

926.11 

926.23 

Ism 

Oi 

929.41 

929.52 

Ism 

Oi 

930.13 

930.26 

Ism 

Hi 

930.61 

930.75 

Ism 

S  VI 

-932 

933.38 

S 

Hi 

936.51 

936.75 

Ism 

Hi 

937.67 

937.80 

Ism 

Svi 

-943 

944.52 

S 

Oi 

948.57 

948.69 

Ism 

Hi 

949.61 

949.74 

Ism 

Oi 

950.76 

950.88 

Ism 

Ni+Oi 

952.26 

Ism 

Ni 

953.30 

953.42 

Ism 

Ni 

953.54 

953.65 

Ism 

Ni 

953.87 

953.97 

Ism 

Ni 

953.97 

954.10 

Ism 

Pii 

962.44 

962.57 

Ism 

Pii 

963.68 

963.80 

Ism 

Ni 

963.87 

963.99 

Ism 

Ni 

964.50 

964.63 

Ism 

Ni 

964.94 

965.04 

Ism 

Oi 

971.61 

971.74 

Ism 

Hi 

972.34 

972.54 

Ism 

Cm 

-976 

977.02 

S 

Cm 

976.88 

977.02 

Ism 

Oi 

976.32 

976.45 

Ism 

Oi 

988.56 

988.65 

Ism 

Si  II 

989.75 

989.87 

Ism 

Nm 

-990 

991.5 

S 

Sim 

997.26 

997.38 

Ism 

Sin 

1020.63 

1020.70 

Ism 

Hi 

1025.37 

1025.72 

S,Ism 

Hi 

1025.68 

1025.72 

S,Ism 

Ovi 

1030.7 

1031.93 

S 

Ovi 

1036.4 

1037.62 

S 

Cii 

1036.24 

1036.34 

Ism 

Oi 

1039.18 

1039.23 

Ism 

An 

1048.15 

1048.22 

Ism 

Siv 

1062.3 

1062.65 

S 

An 

1066.59 

1066.66 

Ism 

Siv 

1072.5 

1072.97 

S 

Nil 

1083.74 

1083.99 

Ism 

Sim 

1109.4 

1108.36 

S 

1109.94 

S 

Sim 

1113.0 

1113.2 

S 

Siv 

1117.3 

1117.76 

S 

Pv 

1117.3 

1117.98 

S 

1127.4 

1128.01 

S 

Si  IV 

1122.0 

1122.5 

S 

Table  2 

(Continued) 


Line 

Identification 

Wavelength 

(A) 

A.0 

(A) 

Origin 

1127.4 

1128.3 

S 

Ni 

1134.08 

1134.16 

Ism 

Ni 

1134.33 

1134.42 

Ism 

Ni 

1134.89 

1134.98 

Ism 

Cm 

1175.0 

1174.9-1176.4 

S 

Notes,  a  =  possibly  with  terrestrial  contamination;  Ism  =  ISM;  S  =  source. 


Table  3 

STIS  Lines 


Line 

Identification 

Wavelength 

(A) 

A.0 

(A) 

Origin 

Cm 

1173.62 

1174.9-1176.4 

S 

Sin 

1190.36 

1190.42 

Ism 

1193.24 

1193.29 

Ism 

Ni 

1199.49 

1199.55 

Ism 

1200.18 

1200.22 

Ism 

1200.67 

1200.71 

Ism 

Sim 

1204.77 

1206.50 

S 

Sim 

1206.45 

1206.50 

Ism 

Hi 

1215.6 

1215.67 

S 

Nv 

1236.9 

1238.82 

S 

1240.3 

1242.80 

S 

Nv 

1238.63 

1238.82 

Ism 

1242.59 

1242.80 

Ism 

Sii 

1250.55 

1250.58 

Ism 

1253.77 

1253.81 

Ism 

1259.45 

1259.52 

Ism 

Sin 

1260.37 

1260.42 

Ism 

Sim 

1297.3 

1294.55 

S 

1296.73 

S 

1298.89 

S 

Oi 

1302.12 

1302.17 

Ism 

Sin 

1304.33 

1304.37 

Ism 

Cii 

1334.47 

1334.53 

Ism 

1334.47 

1334.66 

Ism 

1335.66 

1335.71 

Ism 

Si  IV 

1391.4 

1393.76 

S 

1400.35 

1402.77 

S 

Sin 

1526.62 

1526.71 

Ism 

Civ 

-1545 

1548.20 

S 

1550.78 

S 

Civ 

1547.95 

1548.20 

Ism 

1550.54 

1550.78 

Ism 

Hen 

1637.7 

1640.5 

S 

Alii 

1670.73 

1670.8 

Ism 

Notes.  Ism  =  ISM;  S  =  source;  there  is  possibly  some  atmospheric  absorption 
contamination  for  N  i,  N  ii,  H  i,  and  O  i. 


shows  the  very  good  match  of  the  spectra  in  the  overlap  region, 
while  the  lower  panel  shows  that,  when  normalized,  the  three 
STIS  spectra  track  each  other  with  excellent  fidelity.  The  STIS 
spectra  are  noisy  in  the  overlap  region,  while  the  upper  panel 
of  Eigure  2  shows  that  the  STIS  and  FUSE  spectra  continua 
superpose  accurately.  We  produced  a  combined  spectrum  for 
analysis  by  terminating  the  STIS  spectrum  at  1181.5  A  at 
the  short  wavelength  end,  dividing  it  by  0.98,  and  adding  the 
FUSE  spectrum,  beginning  at  1 1 8 1 .5  A  and  extending  to  shorter 
wavelengths. 

We  suggest,  without  detailed  analysis,  that  the  differences 
among  the  STIS  spectra  could  be  explained  by  slight  variation 
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Figure  2.  Combined  FUSE  and  STIS  spectra.  The  upper  panel  shows  the 
fit  of  the  spectra  in  the  overlap  region  when  small  normalizing  factors  are 
applied  to  the  STIS  spectra.  The  FUSE  spectrum  is  red;  the  first  STIS  spectrum 
(Table  1),  divided  by  0.98,  is  cyan;  the  second  STIS  spectrum,  divided  by  1 . 12,  is 
magenta;  the  third  STIS  spectrum,  divided  by  0.94,  is  brown.  Note  the  accurate 
superposition  of  the  STIS  spectra  with  application  of  normalizing  factors. 

in  the  mass  transfer  rate  from  the  secondary  star.  Whatever  the 
cause  of  the  variation,  the  fact  that  the  normalized  STIS  spectra 
have  the  same  spectral  energy  distribution  (SED)  simplifies 
calculation  of  a  system  model. 

3.  INITIALLY  ADOPTED  SYSTEM  PARAMETERS 

Table  4  (and  notes)  lists  initially  adopted  parameters  and  their 
sources.  Several  of  the  system  parameters  are  poorly  known.  As 
mentioned  earlier,  the  Hipparcos  parallax  is  9.11  ib  1.95  mas, 
which  corresponds  to  a  distance  of  110  zb  30  pc.  Knigge  (2006, 
2007)  determines  M2  =  0.475  Mq  for  P  =  0.20716  days  and 
a  mean  M^d  =  0.7  Mq.  The  mean  M^d  determined  by  Knigge 
is  0.75  Mq,  which  was  inadvertently  truncated  to  0.7  Mq.  This 
change  has  an  imperceptible  effect  on  our  results.  We  adopt  both 
masses  and  note  that  the  resulting  q  =0.68  is  in  agreement  with 
the  observationally  determined  q  ~  0.7  (Hartley  et  al.  2002). 
Panel  et  al.  (2000)  list  a  WD  radius  of  1.091  x  10“^  Rq  for 
a  0.70  Mq  homogeneous  zero  temperature  Hamada-Salpeter 
carbon  model  WD.  In  our  subsequent  study  of  the  observed 
spectra,  we  correct  the  radius  for  the  adopted  WD  Teff;  the 
contribution  of  the  WD  to  the  system  synthetic  spectrum  is 
appreciable.  Tablenotes  to  Table  4  list  alternative  sources  for 
the  adopted  initial  parameters.  These  notes  do  not  include  the 
determination  of  i  ~  25°-30°  (Cowley  &  Crampton  1977). 
This  determination  was  based  on  radial  velocity  observations 
with  a  large  scatter;  the  derived  velocity  amplitude  was  K  = 
74  km  s“\  while  much  more  accurate  data  give  a  value  of 
K  =  165  km  s-i  (Hartley  et  al.  2005). 

Verbunt  (1987)  determined  a  rough  upper  limit  of  E(B  —  V) 
of  0.04  (and  a  preferred  value  of  E(B  —  V)  =  0.0)  for  V3885 
Sgr  based  on  a  study  of  the  2200  A  “bump”  in  International 
Ultraviolet  Explorer  (JUE)  spectra.  Bruch  &  Engel  (1994)  list 
E{B  —  V)  =  0.02,  which  we  adopt. 

A  number  of  studies  have  considered  the  tidal  cutoff  bound¬ 
ary,  rj,  of  accretion  disks  (Paczynski  1977;  Papaloizou  &  Pringle 
1977;  Whitehurst  1988;  Schwarzenberg-Czerny  &  Rozyczka 
1988;  Whitehurst  &  King  1991;  Goodman  1993).  These  au- 


Parameter 

Value 

Parameter 

Value 

Mwd 

0.70^0“ 

i 

65° 

M2 

0.475  Mq'’ 

M 

5.0  X  10-5  y^-l  e 

P 

0.20716071  days ' 

d 

110±30  pcf 

Notes. 

^  Knigge  (2006,  2007)  mean  Mwd;  Hartley  et  al.  (2005),  0.55  Mq  to  0.8  M©; 
Ribeiro  &  Diaz  (2007),  0.3  M©  to  0.9  Mq. 

^  Knigge  (2006,  2007)  calibrated  P  :  M2  relation;  Hartley  et  al.  (2005),  q  =  0.6 
to  0.8;  Ribeiro  &  Diaz  (2007),  M2  =  0.2  M©  to  0.55  Mq. 

^  Ribeiro  &  Diaz  (2007). 

^  Hang  &  Dreschel  (1985),  i  between  60°  and  70°;  Ribeiro  &  Diaz  (2007), 
i  between  45°  and  75°;  Hartley  et  al.  (2005),  i  >  65°. 

^  Linnell  et  al.  (2007)  and  similarity  to  IX  Vel  (Hartley  et  al.  2002). 

^  Hipparcos  value. 

thors  agree  on  ~  0.33D,  where  D  is  the  separation  of  the 
stellar  components.  We  adopt  this  expression  for  the  tidal  cutoff 
radius  of  the  accretion  disk. 

4.  THE  ANALYSIS  PROGRAM:  BINSYN 

Our  analysis  uses  the  program  suite  BINSYN  (Linnell  & 
Hubeny  1996);  recent  papers  (Linnell  et  al.  2007,  2008a) 
describe  its  application  to  CV  systems  in  detail.  Briefiy,  an 
initial  calculation  produces  a  set  of  annulus  models  for  a  given 
WD  mass,  radius,  and  mass  transfer  rate.  This  calculation  uses 
the  program  TLUSTY^  (Hubeny  1988;  Hubeny  &  Lanz  1995). 
The  program  considers  the  radial  and  vertical  structures  of  the 
disk  independently;  the  radial  structure  is  based  on  the  standard 
model  (EKR)  and  so  follows  the  prescribed  relation  between 
local  Teff  and  the  annulus  radius.  The  vertical  structure  is  solved, 
self-consistently,  as  described  by  Hubeny  (1990)  and  Hubeny  & 
Hubeny  (1998).  The  set  of  annulus  models  covers  the  accretion 
disk  from  its  innermost  (WD)  radius  to  just  short  of  its  outermost 
(tidal  cutoff)  radius. 

The  primary  source  of  viscosity  in  CV  accretion  disks  is 
magneto-rotational  instability  (MRI;  Balbus  &  Hawley  1991). 
Hirose  et  al.  (2006)  calculate  magnetohydrodynamics  (MHD) 
models  with  local  dissipation  of  turbulence  and  show  that  the 
vertical  extent  of  an  annulus  is  greater  than  that  in  previous 
models.  Blaes  et  al.  (2006)  show  that  magnetic  support  has 
a  significant  effect  on  synthetic  spectra  of  black  hole  (BH) 
accretion  disk  annuli  and  illustrate  the  effect  in  the  case  of  a  BH 
system  with  Mbh  =  6.62  Mq  and  with  an  adopted  a  =  0.02. 
They  are  able  to  simulate  the  effect  of  magnetic  support  within 
the  TLUSTY  framework  by  adjusting  the  TLUSTY  parameters 
^0  and  and  use  this  simulation  to  calculate  the  effect  on 
synthetic  spectra.  The  primary  effect  is  a  hardening  of  the 
spectrum.  King  et  al.  (2007)  point  out  that  MHD  models  that 
produce  viscosity  via  MRI  require  a  values  that  are  a  factor 
of  10  smaller  than  the  a  =  0. 1-0.4  required  by  observational 
evidence  and  suggest  that  some  caution  is  still  needed  in 
accepting  the  MHD  results.  Our  models  have  used  the  default 
values  ^0  =  ^1  =  0-0,  and  we  assume  that  the  effect  of 
magnetic  support  on  our  annulus  spectra,  with  M^d  =  0.7  Mq, 
is  insignificant. 

Table  5  lists  properties  of  annuli  calculated  with 
TLUSTY(v.202)  for  a  mass  transfer  rate  of  5.0  x  10“^  Mq  yr^. 
This  value  initially  was  chosen  for  test  purposes  because  of  the 
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Table  5 

Properties  of  Accretion  Disk  with  Mass  Transfer  Rate 
M  =  5.0  X  10“^  Mq  yr-i  and  WD  Mass  of  0.70  M© 


^/^wd,0 

Teff 

mo 

logg 

zo 

Ne 

rRoss 

1.36 

53136 

1.221E4 

6.84 

8.27E7 

2.0E17 

2.33E4 

2.00 

47642 

1.491E4 

6.58 

1.43E8 

1.3E17 

2.78E4 

3.00 

38525 

1.456E4 

6.28 

2.43E8 

7.7E16 

3.35E4 

4.00 

32380 

1.356E4 

6.06 

3.49E8 

5.1E16 

3.92E4 

5.00 

28086 

1.259E4 

5.89 

4.60E8 

3.8E16 

4.52E4 

6.00 

24918 

1.176E4 

5.75 

5.76E8 

3.0E16 

5.16E4 

8.00 

20531 

1.043E4 

5.53 

8.14E8 

2.5E16 

6.56E4 

10.00 

17612 

9.439E3 

5.35 

1.06E9 

1.9E16 

8.05E4 

12.00 

15512 

8.664E3 

5.21 

1.32E9 

1.5E16 

9.53E4 

14.00 

13923 

8.040E3 

5.09 

1.58E9 

1.2E16 

1.08E5 

16.00 

12670 

7.526E3 

4.98 

1.85E9 

9.9E15 

1.18E5 

18.00 

11653 

7.092E3 

4.89 

2.12E9 

8.3E15 

1.23E5 

20.00 

10810 

6.721E3 

4.80 

2.38E9 

6.7E15 

1.23E5 

24.00 

9487 

6.115E3 

4.65 

2.92E9 

5.1E15 

1.17E5 

28.00 

8491 

5.638E3 

4.51 

3.57E9 

4.7E15 

1.07E5 

32.00 

7711 

5.251E3 

4.46 

4.90E9 

2.2E14 

9.32E4 

Notes.  Each  line  in  the  table  represents  a  separate  annulus.  A  Shakura  &  Sunyaev 
(1973)  viscosity  parameter  o;  =  0.1  was  used  in  calculating  all  annuli.  The  WD 
radius,  rwd,o?  is  the  radius,  0.01091  R©,  of  a  zero  temperature  Hamada-Salpeter 
carbon  model.  See  the  text  (Section  4)  for  a  discussion  of  the  table  units. 

similarity  of  V3885  Sgr  to  IX  Vel  (Linnell  et  al.  2007);  see 
Section  5  for  a  detailed  discussion.  The  control  file  to  calculate 
a  given  annulus  requires  a  radius  of  the  WD  in  units  of  Rq. 
All  of  the  annuli  used  the  radius  of  a  zero  temperature  WD 
for  a  homogeneous  carbon  Hamada-Salpeter  0.70  Mq  model 
from  Panel  et  al.  (2000).  All  of  the  annuli  are  solar  composi¬ 
tion  H-He  models,  and  the  models  through  r/rwd,o  =  20.00  are 
converged  local  thermodynamic  equilibrium  (LTE)  models.  The 
remaining  annuli  are  LTE-gray  models  (see  the  TLUSTY  Users 
Guide  for  an  explanation).  A  Shakura  &  Sunyaev  (1973)  viscos¬ 
ity  parameter  a  =  0.1  was  used  in  calculating  all  annuli.  Each 
line  in  the  table  represents  a  separate  annulus.  Temperatures  are 
in  Kelvins.  The  column  headed  by  mq  is  the  column  mass,  in 
gm  cm“^,  above  the  central  plane.  The  columns  headed  by  zq 
and  Ne  are,  respectively,  the  height  (cm)  above  the  central  plane 
for  a  Rosseland  optical  depth  of  ~0.7  and  the  electron  density 
(  cm“^)  at  the  same  level.  The  column  headed  by  log  g  is  the  log 
gravity  (cgs  units)  in  the  z-direction  at  a  Rosseland  optical  depth 
of  ~0.7.  The  column  headed  by  trqss  is  the  Rosseland  optical 
depth  at  the  central  plane.  We  pay  attention  to  the  fact  that  the 
annuli  are  optically  thick  to  the  outer  radius  of  the  accretion 
disk. 

Eollowing  the  calculation  of  annulus  models  for  assigned  M^d 
and  M,  program  SYNSPEC^  (v.48;  Hubeny  et  al.  1994)  is  used 
to  produce  a  synthetic  spectrum  for  each  annulus.  We  adopted 
solar  composition  for  all  synthetic  spectra.  The  synthetic  spectra 
include  contributions  from  the  first  30  atomic  species,  assumed 
to  be  in  LTE. 

We  set  SYNSPEC  to  produce  light  intensities  at  10  zenith 
angles  (equally  spaced  in  cos{y)  between  1.0  and  0.1,  where 
y  is  the  angle  between  the  given  direction  and  the  local  nor¬ 
mal)  at  each  calculated  wavelength  for  all  the  individual  an¬ 
nulus  synthetic  spectra.  This  choice  automatically  accounts 
for  wavelength-dependent  limb  darkening  in  calculating  syn¬ 
thetic  spectra  at  different  orbital  inclinations.  BINSYN  calcu¬ 
lates  a  synthetic  spectrum  for  each  of  the  annuli  specified  in 
the  BINSYN  system  model  by  interpolation  among  the  array  of 


^  http://nova.astro.umd.edu 


Table  6 

Temperature  Profile  for  V3885  Sgr  Accretion  Disk  with  the  Mass  Transfer 
Rate  of  M  =  5.0  x  10“^  M©  yr"!  and  WD  Mass  of  0.70  M© 


r/r^d 

Teff 

^/^wd 

Teff 

^/^wd 

Teff 

1.00 

43608 

13.32 

12308 

27.94 

7260 

1.18 

43608 

14.65 

11511 

29.27 

7021 

1.36 

45358 

15.98 

10826 

30.59 

6799 

2.69 

34985 

17.31 

10229 

31.92 

6594 

4.02 

27556 

18.64 

9705 

33.25 

6402 

5.35 

22948 

19.96 

9240 

34.58 

6223 

6.68 

19807 

21.29 

8824 

35.91 

6055 

8.01 

17518 

22.62 

8449 

37.24 

5898 

9.33 

15767 

23.95 

8110 

38.57 

5749 

10.66 

14380 

25.28 

7802 

39.90 

5610 

11.99 

13249 

26.61 

7520 

41.22 

5478 

Notes.  The  WD  radius  rwd  is  the  radius  of  a  60,000  K  0.70  M©  WD  interpolated 
from  Table  4a  of  Panel  et  al.  (2000).  The  Teff  =  43608  K  for  r/rwd  =  1-00 
refers  to  the  inner  edge  of  the  innermost  annulus.  The  same  Teff  for  r/rwd  = 
1.18  refers  to  the  outer  edge  of  the  innermost  annulus  which  coincides  with  the 
inner  edge  of  the  next  annulus.  The  remaining  Teff  values  refer  to  the  inner  edge 
of  the  corresponding  annuli. 

TLUSTY  annulus  models,  with  proper  allowance  for  inclina¬ 
tion. 

BINSYN  models  the  complete  CV  system,  including  the  WD, 
secondary  star,  accretion  disk  face,  and  accretion  disk  rim  as 
separate  entities.  The  model  represents  phase-dependent  and 
inclination-dependent  effects,  including  eclipse  and  irradiation 
effects,  on  all  of  the  system  objects.  The  representation  of  the 
stars  requires  polar  values  and  gravity-darkening  exponents. 
BINSYN  represents  the  accretion  disk  by  a  specified  number  of 
annuli,  where  that  number  typically  is  larger  than  the  number  of 
TLUSTY  annulus  models.  The  accretion  disk  Teff  profile  may 
follow  the  standard  model,  but,  alternatively,  the  profile  may  be 
specified  by  a  separate  input  file. 

BINSYN  has  provision  to  represent  a  bright  spot  on  the 
rim  face,  but  the  low  amplitude  V  light  curve  by  Ribeiro  & 
Diaz  (2007;  Eigure  10)  is  consistent  with  light  variation  of  an 
irradiated  secondary  star,  so  we  do  not  include  a  rim  bright  spot. 
Based  on  the  analysis  by  Smak  (2002),  we  set  the  rim  Teff  equal 
to  the  Teff  of  the  outermost  annulus. 

5.  A  V3885  SGR  SYSTEM  MODEL 

We  adopted  system  parameters  listed  in  Table  4  for  an 
initial  test.  Based  on  the  similarity  to  IX  Vel  (comparable 
periods,  visual  magnitudes,  parallaxes,  CV  class,  spectroscopic 
properties  (Hartley  et  al.  2002)),  we  first  tried  a  standard  model 
accretion  disk  with  the  same  M,  5.0  x  10“^  Mq  yr“^  as  for 
IX  Vel  (Linnell  et  al.  2007),  and,  as  for  IX  Vel,  the  (initial) 
model  included  a  60,000  K  WD.  In  constructing  a  model 
accretion  disk,  it  is  necessary  to  allow  for  the  change  in  WD 
radius  from  the  zero  temperature  model.  We  used  Table  4a 
of  Panel  et  al.  (2000)  to  estimate  the  60,000  K  WD  radius 
at  0.0134  Rq.  Our  BINSYN  model  specified  33  annuli  for  the 
accretion  disk,  with  assigned  radii  and  corresponding  standard 
model  Teff  values  listed  in  Table  6.  The  larger  hot  WD  radius 
than  a  zero  temperature  WD  radius  at  a  given  multiple  of  the 
“WD  radius”  produces  a  lower  Table  6  temperature  than  Table  5 
(e.g.,  at  the  r/r^d  value  6.0  in  Column  1  of  both  tables).  The 
synthetic  spectrum  for  a  given  BINSYN  annulus  is  calculated 
by  interpolation  (temperaturewise)  among  the  Table  5  entries. 
Ideally,  we  could  recalculate  the  whole  series  of  annulus  models 
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(via  TLUSTY)  based  on  the  new  value  of  Rwd,  but  the  change  in 
the  individual  annulus  models  would  be  very  small  and  would 
not  be  warranted  in  view  of  other  uncertainties  (e.g.,  the  system 
parallax). 

Table  5  ends  at  a  Tgff  of  77 1 1  K.  TLUSTY  fails  for  larger  radii, 
even  in  the  gray  model  case.  The  net  effect  is  that  the  calculated 
BINSYN  outer  annuli  have  synthetic  spectra  corresponding  to 
7711  K  even  if  the  specified  standard  model  calls  for  a  lower 
temperature.  The  effect  on  the  synthetic  spectrum  of  failing 
to  follow  the  standard  model  is  negligible.  Adding  a  lower 
temperature  (3500  K)  stellar  synthetic  spectrum  to  the  array  of 
Table  5,  thereby  permitting  interpolation  to  a  lower  temperature 
in  the  outer  BINSYN  annuli,  made  no  detectable  change  in  the 
system  synthetic  spectrum.  As  separate  issues,  the  outer  annuli 
temperatures  should  not  fall  below  6000  K  to  avoid  outbursts 
(Osaki  1996;  Lasota  2001),  and  tidal  dissipation  is  expected  to 
raise  the  outer  accretion  disk  temperature  above  the  standard 
model  value  (Smak  2002). 

The  overall  SED  fit  to  the  combined  spectrum  of  Section  2.3, 
using  initial  parameters,  and  judged  by  visual  estimate,  was  very 
good  with  a  standard  model  accretion  disk  temperature  profile. 
This  was  in  contrast  to  IX  Vel  which  required  modification  of 
the  temperature  profile  from  a  standard  model  to  achieve  an  ac¬ 
ceptable  fit  (Linnell  et  al.  2007).  Assessing  the  quality  of  fit  for 
V3885  Sgr  is  complicated  by  the  presence  of  strong  absorption 
features  of  highly  ionized  material  (Figure  1),  probably  associ¬ 
ated  with  a  wind  (Hartley  et  al.  2002).  At  this  initial  stage,  the 
extreme  synthetic  spectrum  FUV  region  had  too  large  flux:  the 
only  parameter  (either  M  or  Teff^wd)  available  to  adjust  the  fiux 
while  preserving  the  fixed  normalizing  factor  was  the  WD  Tgff . 
(M  could  be  adjusted  but  that  route  produced  detectably  larger 
residuals  longward  of  1500  A.  Small  adjustments  of  M  move 
the  entire  synthetic  spectrum  up  or  down  without  appreciable 
change  in  the  spectral  gradient.  Similarly,  a  change  in  i  does 
not  change  the  spectral  gradient.  We  consider  a  boundary  layer 
in  Section  6.)  We  reduced  the  WD  to  57,000  K  and  found 
this  produced  an  optimum  fit  as  judged  by  eye  estimate.  The 
new  model  neglected  the  very  small  change  in  the  WD  radius 
between  Tgff  =  60,000  K  and  57,000  K,  which  would  have 
affected  the  tabular  radii  in  Table  6. 

Figure  3  shows  the  FUV  fit  to  the  combined  FUSE  and  STIS 
spectra,  dereddened  for  E{B  —  V)  =  0.02.  We  pay  attention 
to  the  broad  absorption  lines  which  we  attribute  to  a  wind/ 
chromosphere  and  that  affect  a  major  part  of  the  FUSE  spectrum. 
Note  the  deep  ISM  Lyman  series  lines,  the  broad  and  deep 
absorption  by  O  i  and  N  iii  near  990  A,  the  deep  and  possible 
multi-component  line  near  1025  A,  and  the  broad  lines 
of  Sim,  Pv,  Siv,  and  Siiv  in  the  1110-1140  A  interval. 
The  overlapping  and  possibly  multicomponent  lines  depress 
the  apparent  continuum  over  the  entire  interval  from  970  A  to 
1170  A. 

The  normalizing  factor  to  divide  the  synthetic  spectrum  to 
superpose  it  on  the  observed,  dereddened  combined  FUSE  and 
STIS  spectrum,  corresponding  to  an  exact  distance  of  110  pc, 
is  1.1473  X  10"^^  cm^.  We  stress  that  this  is  not  an  adjustable 
parameter.  In  Section  5.1,  we  consider  the  effects  of  varying 
the  “exact”  distance  corresponding  to  dilcr  in  the  Hipparcos 
parallax.  The  spectrum  at  the  bottom  is  the  contribution  of  a 
57,000  K  WD.  For  simplicity,  we  assume  that  the  WD  is  rotating 
synchronously;  we  do  not  fit  any  observed  absorption  lines  in 
detail.  The  program  assumes  that  the  accretion  disk  blocks  the 
view  of  the  lower  half  of  the  WD.  In  addition,  BINSYN  allows 
for  the  partial  visibility  of  that  part  of  the  WD  not  obscured  by 
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Figure  3.  FUV  fit  of  system  synthetic  spectrum  to  FUSE  and  STIS  combined 
spectrum  (gray  plot)  of  V3885Sgr.  Note  the  extensive  absorption  by  local  highly 
ionized  material.  The  lower  spectrum  is  the  contribution  of  a  57,000  K  WD;  the 
top  synthetic  spectrum  is  the  system  synthetic  spectrum  including  the  WD  and 
a  M  =  5.0  X  10“^  Mq  yr“^  standard  model  accretion  disk.  See  the  text  for  a 
discussion. 

the  accretion  disk  and  above  the  observer’s  horizon  and  includes 
wavelength-dependent  limb  darkening.  The  spectral  resolution 
of  all  the  synthetic  spectra  is  0.02  A  (Section  5.2.).  The  adopted 
orbital  phase  for  the  synthetic  spectra  was  0.75,  in  agreement 
with  the  observed  spectra.  Section  2.  Adopting  an  orbital  phase 
of  0.25  produced  a  clear  misalignment  of  the  calculated  960  A 
peak  with  the  observed  one. 

We  check  the  apparent  best  fit,  by  eye  estimate,  with  a  reduced 
chi-squared  (x^)  test.  We  did  not  smooth  the  observed  spectrum 
for  this  test.  The  presence  of  unmodeled  wind/chromosphere 
features  in  the  observations  produces  large  residuals  (especially 
in  the  FUV,  Figure  3).  A  x^  test  of  our  model  fit  to  the  observed 
continuum  requires  masking  those  spectral  regions  excluded 
from  the  fit  (i.e.,  the  wind/chromosphere  features).  We  initially 
masked  spectral  lines  by  hand,  including  the  deep  Lyman  lines. 
(The  Lyman  lines  in  the  model,  mostly  from  the  accretion 
disk,  are  washed  out  by  Keplerian  rotation.)  This  first  stage 
masking  produced  a  major  reduction  (from  154.39  to  7.38)  in 
X^  as  compared  with  the  unmasked  spectrum.  There  remained 
a  major  “valley”  in  a  residuals  plot,  described  in  the  following 
paragraph,  with  the  largest  (negative)  residuals  near  1 100  A.  Our 
previous  discussion  has  identified  this  effect  with  overlapping 
absorption  features  from  a  wind/chromosphere  which  depresses 
the  entire  spectral  interval  over  the  approximate  range  970  A- 
1170  A.  BINSYN  has  no  present  capability  of  modeling  a 
wind/ chromosphere,  and  further  progress  requires  masking  the 
affected  spectral  interval.  We  therefore  masked  the  entire  region 
from  970  A  to  1170  A,  producing  a  further  large  reduction  in 
X^  (from  7.38  to  2.65). 

Figure  4  shows  plots  of  the  residuals  in  the  three  cases. 
The  top  plot,  displaced  upward  by  12.0  ordinate  units,  is  the 
residuals  for  the  unmasked  case.  The  middle  plot,  displaced 
upward  by  4.0  ordinate  units,  is  the  partially  masked  case.  The 
lower  (undisplaced)  plot  is  the  final  set  of  residuals  adopted  for 
this  study.  In  each  case  a  horizontal  line  marks  the  undisplaced 
reference.  The  change  from  the  top  to  middle  plot  mostly  results 
from  masking  of  individual  spectral  features,  such  as  the  Lyman 
lines,  the  Civ  kl548  feature,  etc.  The  STIS  spectrum  appears 
to  have  a  defect  beginning  at  1710  A,  and  the  three  STIS 
spectra  show  more  variability  longward  of  the  Hen  kl640  line 
than  in  other  spectral  regions.  We  terminated  our  calculated  fit 
just  longward  of  the  Hen  kl640  line.  The  middle  plot  clearly 
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Table  7 

Values  of  x  ^  as  Function  of  Model  Parameters 


M  Mq  yd) 

7wd(K) 

/(deg.) 

M  (10-9  Mq  yd) 

7wd(K) 

/(deg.) 

x^ 

4.0 

70,000 

65 

7.51 

5.0 

57,000 

64 

3.95 

4.0 

75,000 

65 

5.60 

5.0 

57,000 

66 

4.43 

4.5 

57,000 

65 

6.86 

5.0 

57,000 

70 

39.73 

4.5 

60,000 

65 

5.05 

5.2 

40,000 

65 

4.69 

4.5 

65,000 

65 

3.93 

5.2 

45,000 

65 

3.14 

4.5 

70,000 

65 

3.14 

5.2 

50,000 

65 

2.91 

4.5 

75,000 

65 

3.97 

5.2 

52,000 

65 

2.73 

4.8 

52,000 

65 

5.22 

5.2 

57,000 

65 

3.04 

4.8 

57,000 

65 

3.42 

5.2 

60,000 

65 

4.41 

4.8 

60,000 

65 

3.00 

5.5 

35,000 

65 

4.34 

4.8 

65,000 

65 

3.40 

5.5 

40,000 

65 

3.08 

4.8 

70,000 

65 

4.26 

5.5 

45,000 

65 

3.19 

5.0 

45,000 

65 

4.54 

5.5 

50,000 

65 

3.24 

5.0 

50,000 

65 

4.25 

5.5 

52,000 

65 

3.93 

5.0 

52,000 

65 

3.39 

5.5 

57,000 

65 

5.93 

5.0 

57,000 

65 

2.65 

6.0 

35,000 

65 

4.03 

5.0 

60,000 

65 

3.44 

6.0 

40,000 

65 

6.48 

5.0 

65,000 

65 

5.80 

6.0 

45,000 

65 

4.13 

5.0 

57,000 

60 

49.94 

6.0 

50,000 

65 

9.88 
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Figure  4.  Plots  of  residuals  from  the  visually  selected  best  model.  The  top  plot, 
displaced  upward  by  12.0  ordinate  units,  is  the  residuals  in  the  unmasked  case. 
The  middle  plot,  displaced  upward  by  4.0  ordinate  units,  is  the  residuals  for  a 
partially  masked  observed  spectrum.  The  bottom  plot  is  the  residuals  for  the 
hnally  adopted  masked  case.  Note  the  three  horizontal  lines  marking  residual 
values  of  0.0  for  the  three  cases.  See  the  text  for  a  discussion. 


shows  the  residual  depression  after  individual  spectral  features 
have  been  masked,  and  the  bottom  plot  shows  the  overall 
improvement  in  residuals  after  masking  the  970-1 170  A  region. 
The  improvement  is  even  more  apparent  on  individual  plots  with 
extended  ordinate  axes. 

The  lower  plot  still  has  fairly  large  outliers.  Although  further 
selective  masking  could  be  done,  we  consider  it  dangerous 
to  go  very  far  in  suppressing  observations  that  disagree  with 
the  model.  Residual  absorption  features  could  be  masked  but 
there  are  positive  outliers  that  then  would  bias  the  solution. 
We  anticipate  that  it  will  be  possible  to  locate  an  optimum 
model  with  the  chosen  masked  set  of  observations.  The  final 
masked  version  of  the  observed  combined  spectra  includes  4332 
tabulated  wavelengths  and  corresponding  flux  values. 

It  is  particularly  unfortunate  for  this  analysis  that  most  of 
the  FUSE  spectrum  must  be  masked,  but  it  is  unavoidable:  the 
masked  region  clearly  is  strongly  affected  by  the  broad  and 
deep  absorption  lines  of  high  excitation  species  as  described 


in  Section  2.3.  The  high  excitation  involved  excludes  the  lines 
from  production  by  the  accretion  disk  itself  or  the  WD.  Eurther, 
ignoring  the  high  excitation  aspect  for  the  moment,  it  would 
not  be  possible  to  simulate  the  depressed  970-1170  A  spectral 
interval  by  adopting  a  modified  accretion  disk  Tgff  profile.  Our 
experience  with  non-standard  model  accretion  disks  is  that  the 
calculated  SED  profile  can  have  its  slope  changed  over  an 
extended  spectral  region,  but  it  is  not  possible  to  put  a  “hole”  in 
the  SED. 

We  calculated  values  for  a  grid  of  34  models  in  the  M, 
T'eff.wd  plane,  centered  on  the  best  visual  estimate  model  fit  to 
the  masked  observations.  Table  7  lists  those  values  plus 
values  for  four  additional  i  values  (to  be  discussed  separately). 
The  aim  is  to  provide  a  sufficiently  extended  grid  that  a  subset 
of  x^  values  defines  a  x^  contour  that  encloses  the  region  of 
minimum  x^,  given  the  adopted  values  of  other  parameters.  The 
calculation  of  a  given  model  is  a  time-consuming  process,  and 
a  compromise  is  necessary  between  total  computational  time 
and  density  of  grid  points.  Our  computed  grid  is  sparse,  but  our 
judgment  is  that  it  meets  our  needs,  given  the  uncertainties  of 
separately  adopted  parameters.  We  have  used  the  IDE  routine 
CONTOUR  to  produce  contour  plots  for  the  i  =  65°  values  of 
Table  7.  The  M  =  5.0  x  IQ-^  Mq  yr-\  reff,wd  =  57,000  K 
model  has  the  smallest  tabular  x  ^  value,  and  we  take  that  model 
to  be  our  solution  for  V3885  Sgr.  We  discuss  this  choice  and  the 
topic  of  grid  density,  including  its  adequacy,  in  more  detail  in 
Section  6. 

We  have  sampled  the  M,  Teff^wd  plane  for  the  case  M^d  = 
0.8  Mq  with  two  new  models:  M  =  5.0  x  10“^MQyr“^ 
and  Teff^wd  =  57,000  K,  producing  x^  =  3.06,  and  M  = 
5.0  X  10“^  Mq  yr~^  and  Teff^wd  =  52,000  K,  producing  x^  = 
3.81.  The  corresponding  x^  values  from  Table  7  are  2.65  and 
3.39,  respectively.  A  preferred  value  of  M^d  is  not  possible  by 
X  ^  analysis  of  the  existing  data.  The  SED  accretion  disk  model 
is  not  a  strong  function  of  M^d  for  the  change  from  0.7  Mq 
to  0.8  Mq.  Said  differently,  Eigure  5  remains  roughly  constant, 
with  a  small  systematic  shift  for  changes  in  M^d  of  about  0.1  or 
less. 

An  improved  value  of  M^d  may  be  possible  with  further 
optical  region  spectroscopy.  Ribeiro  &  Diaz  (2007)  find  that 
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Figure  5.  Contour  plot  showing  values  of  in  the  M,  Teff^wd  plane  for  the 
system  models  tabulated  in  Table  7.  An  asterisk  marks  the  value  for  the 
preferred  solution.  See  the  text  for  a  discussion. 

the  Ha  emission  line  arises  from  irradiation  of  the  secondary, 
while  the  (doubled)  A.6678  Hei  emission  line  associates  with 
the  accretion  disk.  If  an  accretion  disk  model  that  produces 
secondary  component  emission  lines  were  available,  a  fit  to  the 
He  I  line  could  constrain  i. 

Table  7  includes  four  models  for  different  assumed  values  of 
i  and  with  other  parameters  the  same  as  the  preferred  model. 
There  is  a  very  sharp  minimum  of  at  the  adopted  i  =  65°. 
At  this  inclination,  there  is  an  eclipse  of  the  outer  accretion 
disk  at  orbital  phase  0.0  with  photometric  effects  much  smaller 
than  the  predicted  variation  due  to  irradiation  of  the  secondary. 
Increasing  the  inclination  to  70°  greatly  deepens  the  depth  of 
accretion  disk  eclipse  and  barely  fails  to  eclipse  the  WD;  it 
would  produce  light-curve  effects  that  are  not  observed  (Ribeiro 
&  Diaz  2007,  Figure  10).  However  the  observed  light  variation 
is  complex  and  may  include  a  contribution  from  a  rim  bright 
spot.  An  inclination  of  71°  would  eclipse  the  WD  and  produce 
a  major  light-curve  variation.  It  is  not  possible  to  use  the 
i  -dependent  variation  to  determine  i.  The  calculated  system  fiux 
is  a  strong  function  of  i  because  of  the  changing  projection 
of  the  flux-dominating  accretion  disk  on  the  plane  of  the  sky. 
The  fit  of  the  model  spectrum  to  the  observations  requires  a 
normalizing  factor  that  is  given  by  the  Hipparcos  parallax. 
The  sharp  minimum  of  x^  for  i  =  65°  results  from  the  joint 
consistency  of  that  parallax  and  the  adopted  i  value.  Actual 
determination  of  i  would  require  observational  effects  intrinsic 
to  the  system,  i.e.,  eclipses. 

Until  a  model  that  includes  simulation  of  the  wind/ 
chromosphere  is  available,  it  is  doubtful  that  substantial  fur¬ 
ther  reduction  in  x^  can  be  achieved.  Based  only  on  the  x^ 
plot.  Figure  5,  we  estimate  the  parameter  errors  in  M,  Teff^wd. 
and  i  as  2.0  x  10“^^  Mq  yr“^  3000  K,  and,  from  Table  7,  2?0, 
respectively.  However,  shifts  in  M  to  match  the  uncertainty  in 
the  Hipparcos  parallax  (Section  5.1.)  require  a  larger  M  error 
of  2.0  X  10“^  Mq  yr“^  which  we  adopt  for  the  system  solution. 
Further,  in  recognition  of  the  sparse  grid  of  models  and  their 
associated  x^  values,  we  increase  the  estimated  wd  error  to 
5000  K. 

Our  final  system  parameters  are  in  Table  8.  Where  appropriate 
we  have  provided  estimated  error  limits.  Figure  6  shows  the 
final  model  compared  with  the  combined  observations  (not 
the  masked  observations).  The  blue  line  is  the  contribution  of 


Parameter 

Value 

Parameter 

Value 

Afwd 

0.70  ±0.1  Mo 

7;ff,s(pole) 

3650  ±  100  K 

Ms 

0.475  Mq 

reff,s(point) 

2214  K 

M 

5.0  ±2.0  X  10-^  Mo  yr-i 

7;ff,s(side) 

3598  K 

P 

0.2071607  days 

Teff,s(back) 

3489  K 

D 

1.55396  Rq 

rXpole) 

0.209  Rq 

Qwd 

116.0 

(point) 

0.296  Rq 

3.2053112 

reside) 

0.218  Rq 

i 

65°  ±  2° 

r^back) 

0.239  Rq 

^eff,wd 

57,000  ±  5000  K 

log  g,(pole) 

4.13 

^ wd 

0.0134  Rq 

log  g,  (point) 

-0.390 

log  gwd 

8.02 

log  gXside) 

4.65 

^wd 

1.0 

log  g,(back) 

4.48 

A. 

0.6 

ra 

0.56  Rq 

^wd 

0.25 

rb 

0.01091  Rq 

hs 

0.08 

H 

0.10  Rq 

Notes,  “wd”  refers  to  the  WD;  5  refers  to  the  secondary  star.  D  is  the  component 
separation  of  centers,  Q  is  a  Roche  potential.  Temperatures  are  polar  values,  A 
values  are  bolometric  albedos,  and  ^  values  are  gravity-darkening  exponents.  Va 
specifies  the  outer  radius  of  the  accretion  disk,  set  at  the  tidal  cutoff  radius,  and 
is  the  accretion  disk  inner  radius,  as  determined  in  the  final  system  model, 
while  H  is  the  semi-height  of  the  accretion  disk  rim. 


Figure  6.  As  in  Figure  3  but  for  the  full  wavelength  range  of  the  combined 
spectra.  The  synthetic  spectrum  of  the  accretion  disk  alone  is  the  next-to-top 
synthetic  spectrum.  The  uppermost  system  synthetic  spectrum  is  the  sum  of  the 
accretion  disk  and  the  WD  (blue  line).  The  contributions  of  the  secondary  star 
and  accretion  disk  rim  are  negligible. 

the  57,000  K  WD.  The  lower  of  the  two  black  lines  is  the 
contribution  of  the  accretion  disk,  and  the  upper  black  line  is 
the  system  model. 

5.7.  The  Effects  of  Hipparcos  Parallax  Uncertainty 

It  is  possible  to  test  the  effect  of  uncertainty  in  the  Hipparcos 
parallax  by  changing  the  normalizing  divisor  for  superposing 
the  synthetic  spectrum  on  the  combined  FUSE  and  STIS 
spectra.  Changing  the  parallax  by  \a  to  7.16  mas  changes 
the  normalizing  divisor  from  1.147  x  10"^^  cm^  to  1.857  x 
10"^^  cm^.  The  synthetic  spectrum  is  displaced  downward  by 
approximately  two  fiux  units  in  Figure  3  (the  1000  A  peak 
moves  from  an  ordinate  value  of  10.  to  8).  Compensating 
for  the  change  would  require  increasing  M  to  perhaps  6.5  or 
7.0  X  10“^  Mq  yr“^  Correspondingly,  changing  the  parallax  to 
1 1.06  mas  changes  the  normalizing  divisor  to  7.784  x  10"^^  cm^ 
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Parameter 

Value 

N(H) 

1.0  X  10^^  cm-2 

DjH 

2  X  10-5 

N(H2) 

1  X  10^^  cm-^ 

Vel. 

—50.0  km  8“^ 

Temp. 

250  K 

Metallicity 

1.0  (1.0  =  solar) 

and  displaces  the  synthetic  spectrum  upward  by  about  two  flux 
units,  from  a  1000  A  peak  of  10  to  12.  Compensating  would 
require  decreasing  M  to  perhaps  4.0  or  4.5  X  10“®  Mo  yr-'. 

5.2.  Correction  for  the  ISM 

Recent  papers  have  included  corrections  for  line  absorption 
by  the  ISM  (Godon  et  al.  2008;  Linnell  et  al.  2008a).  A  custom 
spectral  fitting  package  is  used  to  estimate  the  temperature  and 
density  of  the  interstellar  absorption  lines  of  atomic  and  molec¬ 
ular  hydrogen.  The  ISM  model  assumes  that  the  temperature, 
bulk  velocity,  and  turbulent  velocity  of  the  medium  are  the  same 
for  all  atomic  and  molecular  species,  whereas  the  densities  of 
atomic  and  molecular  hydrogen,  and  the  ratios  of  deuterium  to 
hydrogen  and  metals  (including  helium)  to  hydrogen  are  ad¬ 
justable  parameters.  The  model  uses  atomic  data  of  Morton 
(2000,  2003)  and  molecular  data  of  Abgrall  (2000).  The  molec¬ 
ular  hydrogen  transmission  values  have  been  checked  against 
those  of  McCandliss  (2003).  The  calculation  requires  fairly  high 
spectral  resolution,  and  we  have  calculated  all  of  the  TLUSTY- 
SYNSPEC  synthetic  spectra  at  a  resolution  of  0.02  A;  we  used 
the  same  resolution  in  calculating  the  system  synthetic  spectra. 
The  ISM  parameters  are  in  Table  9.  Eigure  7  shows  the  cor¬ 
rected  version  of  part  of  Eigure  6.  The  ISM  model  produces  a 
good  representation  of  the  observed  deep  cores  of  the  Lyman 
series  lines  and  the  Lyman  cutoff,  as  well  as  some  lines  in  other 
spectral  regions,  and  clearly  shows  that  the  ISM  produces  the 
observed  Lyman  lines.  The  broad,  high  excitation  absorption 
regions  from  the  accretion  disk  corona,  chromosphere,  or  wind, 
identified  in  Tables  2  and  3,  remain  as  prominent  residuals.  As 
stated  above,  the  Lyman  lines  were  masked  in  the  fit  to  the 
model,  and  the  spectral  region  shortward  of  the  Lyman  limit 
was  masked  so  only  the  Lyman  region  continuum  participated 
in  the  model  solution.  This  was  proper  since  a  fit  to  the  ISM  does 
not  affect  the  choice  of  the  model  parameters  listed  in  Table  8. 
There  are  small  positive  residuals  from  the  model  near  the  short 
wavelength  end.  This  could  be  interpreted  as  evidence  for  an 
unmodeled  very  high  temperature  contribution,  discussed  in  the 
following  section. 

6.  DISCUSSION 

In  a  well-known  paper,  Wade  (1988)  showed  that  neither 
steady  state  model  accretion  disks  based  on  Planck  functions  nor 
stellar  model  atmospheres  could  simultaneously  fit  the  colors 
and  absolute  luminosities  of  a  set  of  NL  systems  (V3885  Sgr 
was  included  in  the  study).  Tomographic  analyses  of  accretion 
disks  (Rutten  et  al.  1992)  show  clear  departures  from  standard 
model  temperature  profiles  (LKR).  Recent  analyses  of  the  NL 
systems  MV  Lyr,  IX  Vel,  QU  Car,  and  UX  UMa  (Linnell  et  al. 
2005,  2007, 2008a,  2008b)  showed,  in  each  case,  that  a  standard 
model  accretion  disk  synthetic  spectrum  could  not  accurately  fit 
observed  spectra.  These  results  suggest  that  physical  conditions 
in  possibly  all  NL  systems  are  so  complex,  with  winds,  coronae. 


Figure  7.  Model  of  Figure  6  but  with  the  ISM  correction.  Note  the  cutoff  at 
912  A. 


iron  curtains,  etc.,  that  an  accurate  fit  by  a  standard  model 
synthetic  spectrum  is  too  much  to  expect.  Consequently,  the 
most  interesting  result  of  this  study  is  that  a  standard  model 
synthetic  spectrum  plus  WD  accurately  fit  the  observations 
(X^  =  2.65)  after  masking  the  effects  of  a  wind/chromosphere. 
This  conclusion  is  particularly  true  in  light  of  the  fact  that  the 
system  distance  is  known  so  the  fit  is  in  absolute  flux  units.  It 
is  worth  emphasizing  that  the  model  presented  here  includes  a 
measurable  contribution  from  a  hot  source,  in  addition  to  the 
accretion  disk,  sufficiently  large  that  limits  can  be  set  on  the 
hot  source  Teff.  Assuming  that  the  hot  source  is  a  WD  (see 
below),  this  is  the  first  time  the  effective  temperature  of  a  hot 
WD  has  been  constrained  during  the  high  brightness  state  of  a 
NL  variable. 

As  discussed  by  LKR,  half  of  the  potential  energy  liberated  in 
the  fall  of  the  mass  transfer  stream  from  the  LI  point  (essentially 
from  infinity)  to  the  WD  surface  appears  as  radiated  energy. 
Energy  conservation  requires  an  accounting  for  the  other  half, 
and  the  presence  of  a  hot  boundary  layer  (BL)  is  a  common 
prescription.  The  problem  is  that  BLs  with  the  prescribed 
properties  are  not  observed  (Cordova  1995).  The  absence  of 
an  observed  but  predicted  (LKR)  BL  in  high  M  CV  systems  has 
an  extensive  history:  (Lerland  1982;  Kallman  &  Jensen  1985; 
Patterson  &  Raymond  1985;  Hoare  &  Drew  1991, 1993;  Vrtilek 
et  al.  1994;  Idan  &  Shaviv  1996).  One  proposed  explanation 
is  that  the  WD  is  rotating  close  to  Keplerian  rotation,  with 
no  BL  predicted.  However,  Hubble  Space  Telescope  (HST) 
measurements  starting  with  Sion  et  al.  (1994)  have  measured 
rotational  WD  velocities  too  low  to  explain  “missing”  BLs.  An 
alternative  explanatory  theme  is  that  the  absence  of  an  observed 
BL  associates  with  a  wind,  and  in  V3885  Sgr  a  wind  is  clearly 
present  (Hartley  et  al.  2002,  2005). 

BINSYN  has  an  option  to  simulate  a  BL  by  specifying  the 
innermost  annulus  radial  thickness  and  Tgff.  We  applied  the 
theory  of  Godon  et  al.  (1995)  to  obtain  those  two  parameters. 
This  theory  derives  a  BL  with  a  Teff  ~  100,000  K  and  radial 
thickness  of  about  0.1  This  BL  Tgff  is  consistent  with 
a  value  from  wind  ionization  (Hoare  &  Drew  1991)  and  the 
theory  of  Idan  &  Shaviv  (1996).  The  Godon  et  al.  (1995)  model 
adopted  M^d  =  1.0  Mq  and  7kd,eff  =  20,000  K.  Addition  of  the 
BL  made  a  nearly  negligible  change  in  our  system  SED.  There 
was  a  barely  detectable  increase  in  the  system  flux  at  the  short 
wavelength  limit  over  the  SED  of  our  original  model,  producing 
a  slightly  better  fit  in  that  spectral  region.  The  evidence  in 
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favor  of  a  BL  is  inconclusive;  a  model  including  a  100,000  K 
BL  and  a  57,000  K  WD  is  consistent  with  the  observations. 
Popham  &  Narayan  (1995)  calculate  BL  models  and  illustrate 
their  theory  for  the  case  M  =  10“^  Mq  yr“^  and  zero  WD 
rotation,  leading  to  a  somewhat  higher  BL  Jeff.  In  a  discussion 
among  the  co-authors,  Godon  noted  that  the  BL  Teff  might  be  as 
high  as  200,000  K  but  probably  would  not  be  visible  at  i  =65° 
because  of  an  inflated  inner  disk — an  additional  reason  for  not 
including  a  BL.  A  model  including  the  100,000  K  BL  but  with 
our  WD  contribution  suppressed  (i.e.,  requiring  the  hot  source 
to  be  the  BL  exclusively)  would  be  strongly  inconsistent  with 
the  observations.  We  are  unaware  of  any  BL  models  that  treat 
MHD  effects. 

Is  it  possible  to  exclude  a  different  source  for  the  WD  contri¬ 
bution  of  our  present  model?  We  seek  some  spectral  feature  that 
could  distinguish  the  WD  contribution  from  a  different  prospec¬ 
tive  source.  In  Figure  6,  the  WD  Ly  y  absorption  line  appears 
to  make  an  important  contribution  to  the  system  synthetic  spec¬ 
trum.  A  different  smoothly  varying  source,  such  as  a  black  body 
(BB)  or  power-law  source  might  produce  a  much  poorer  fit  to 
the  observed  Ly  y  wings  and  so  be  excludable.  We  tested  this 
possibility  by  replacing  the  WD  synthetic  spectrum  by  a  57,000 
K  BB  with  the  same  size  as  the  replaced  WD.  The  existing  model 
uses  light  intensities  at  10  zenith  angles  (Section  4),  produced 
by  SYNSPEC,  for  each  annulus.  Use  of  the  BB  spectrum  re¬ 
quired  a  change  to  flux-based  synthetic  spectra  with  an  adopted 
wavelength-independent  limb  darkening  coefficient;  we  chose 
0.6.  The  system’s  synthetic  spectrum’s  SED  was  essentially  un¬ 
changed  but  required  a  slight  adjustment  of  the  normalizing 
factor  to  reproduce  the  same  fit  as  with  the  original  model, 
likely  a  result  of  the  change  to  a  wavelength-independent  limb 
darkening  coefficient.  The  fit  at  the  short  wavelength  limit  was 
slightly  better,  but  it  is  uncertain  whether  this  arises  from  the 
BB  replacement  or  the  change  in  limb  darkening.  In  any  event, 
the  Ly  y  line  was  slightly  shallower  in  the  new  simulation,  ar¬ 
guing  weakly  in  favor  of  the  presence  of  a  WD,  but  the  fit  to  the 
wings  was  essentially  the  same.  The  proposed  spectral  feature 
test  fails,  and  we  are  unable  to  exclude  a  BB  as  the  source  of 
the  required  hot  addition  to  the  accretion  disk.  The  SED  of  a 
power-law  source  would  differ  from  that  of  a  BB  source  and 
should  produce  a  poorer  lit  to  the  observed  spectrum;  without 
explicit  test  we  discount  this  type  of  the  source  as  a  reasonable 
candidate. 

To  answer  the  question  “What  is  the  source  of  the  hot  addition 
to  the  accretion  disk?”  the  most  reasonable  working  hypothesis 
is  that  the  hot  source  is  a  WD.  A  WD  is  known  to  be  present  (it  is 
a  CV  system)  and  the  inclination,  /,  is  such  that  the  WD  should 
be  visible.  A  spectral  source  substitute  for  a  WD  (e.g.,  BB) 
would  need  not  only  to  explain  why  the  WD  spectrum  is  not 
seen  but  also  why  a  BB  source  of  the  same  Tgff,  size  (which 
sets  the  radii  of  the  accretion  disk  annuli  and  thus  their  Teff 
values  (Table  6)),  and  physical  location  (exactly  centered  on  the 
accretion  disk)  is  a  reasonable  substitute. 

The  adequacy  of  our  set  of  34  models  populating  the  x^(M, 
7iff,wd)  plane  to  determine  an  optimum  model  requires  justifi¬ 
cation.  At  issue  is  the  question  “How  dense  a  grid  is  adequate?” 
Phrased  differently,  “As  we  produce  a  denser  and  denser  grid, 
and  keeping  in  mind  the  commitment  of  computer  time,  how  do 
we  know  when  to  stop?”  A  very  important  point  is  that  our  so¬ 
lution  adopts  values  for  the  parameters  M(wd),  M(sec),  period, 
/,  parallax,  and  E{B  —  V).  A  solution  with  a  different  value  of 
any  one  of  these  would  require  recalculation  of  an  entire  new 
plane  of  x^(M,  Teff^wd)  values,  like  Eigure  5,  and  search  of  it  for 


the  minimum  x^.  This  process  would  be  required,  for  example, 
if  we  study  the  implications  of  changing  the  adopted  parallax 
by  one  standard  error.  The  solution  we  have  found  from  our  set 
of  x^  values  is  valid  only  under  the  assumption  that  the  other 
fixed  parameter  values  are  accurate. 

Our  Table  7  values  for  a  specific  M  and  variable  Teff^wd,  or 
specific  Teff^wd  and  variable  M,  extend  to  a  x^,  when  rounded  to 
the  nearest  half-integer,  of  at  least  4.0.  We  consider  this  value  to 
be  enough  larger  than  the  minimum  of  2.65,  and  the  resolution 
in  M  and  Teff^wd  sufficiently  fine  to  provide  reasonable  assurance 
that  the  (ideal)  global  minimum  has  been  encircled.  Eigure  5, 
produced  with  the  IDL  routine  CONTOUR,  presents  contours  of 
constant  x^  including  closed  contours  that  encircle  the  tabular 
minimum. 

We  reiterate  the  point  that  possibly  our  most  interesting  result 
is  the  accurate  fit  to  a  standard  model.  We  note  that  a  visual 
estimate  (“chi  by  eye”)  selected  the  best  model. 

Hartley  et  al.  (2002)  discuss  the  wind  in  V3885  Sgr  and  point 
to  the  lack  of  correlation  between  wind  activity  and  system 
luminosity,  an  inconsistency  with  expectation  for  line-driven 
winds.  The  theory  of  line-driven  winds  has  been  applied  to  CV 
systems  by  Pereyra  et  al.  (2006)  and  Proga  (2003a,  2003b).  See 
the  review  by  Drew  &  Proga  (2000).  Proga  (2003a)  considers 
both  radiation-driven  and  MHD  effects  and  shows  that  hybrid 
niodels  that  include  both  MHD  and  line-driving  produce  larger 
4^wind  values,  of  order  10“^^  Mq  yr“^  but  Proga  (2003b)  shows 
that  a  pure  line-driven  wind  produces  an  improved  fit  to  observed 
line  shapes  but  with  a  lower  M^ind-  As  Hartley  et  al.  (2002)  point 
out,  the  lack  of  correlation  between  wind  activity  and  system 
luminosity  may  argue  in  favor  of  MHD  effects.  In  either  case, 
the  theoretically  derived  values  of  M^ind  are  far  smaller  than  the 
mass  transfer  M,  and  essentially  all  of  the  mass  transfer  stream 
must  be  deposited  on  the  WD. 

Using  Equation  (2)  of  Sion  (1995),  the  Teff  of  a  compression- 
ally  heated  WD  is  given  by  Teff  =  [(/GMwdM)/(47rR^^a)]^-^^, 
where  /  is  the  fraction  of  the  accretion  energy  applied  to  com- 
pressional  heating  and  is  in  the  range  0.15-0.20.  Eor  /  =  0.15 
and  the  parameters  of  Table  8,  we  calculate  Teff  =  52,500  K, 
while  for  /  =  0.20,  we  And  T^ff  =  56,414  K,  well  within  the 
estimated  error  limit  of  our  57,000  K  WD. 

This  study  raises  an  obvious  question:  why  does  a  standard 
model  accretion  disk  represent  V3  885  Sgr  accurately  while  other 
apparently  comparable  systems  violate  the  standard  model? 
To  rephrase  the  question:  what  equation  or  equations  defining 
the  standard  model  do  the  nonconforming  systems  violate? 
The  standard  model  equations  (e.g.,  EKR  Equation  (5.68)) 
adopt  axial  symmetry  and  hydrostatic  equilibrium,  avoid  MHD 
considerations  (the  MRI  (Balbus  &  Hawley  1991)  source  of 
viscosity),  and  assume  a  constant  a  for  the  entire  accretion  disk, 
corresponding  to  a  constant  local  M  within  the  accretion  disk. 
The  observed  presence  of  spiral  waves  in  V3885  Sgr  (Hartley 
et  al.  2005)  violates  the  axial  symmetry  condition  even  in  a 
system  that  fits  the  standard  model.  The  impact  of  the  mass 
transfer  stream  produces  a  departure  from  axial  symmetry. 
These  phenomena  illustrate  the  point  that  the  standard  model 
is  a  simplified  but  very  useful  idealization.  Just  as  weather  in 
the  Earth’s  atmosphere,  from  physically  complex  phenomena, 
produces  transient  but  extended  local  departures  from  the 
ideally  accurate  hydrostatic  equilibrium,  we  speculate  that 
transients  of  unknown  but  physically  complex  origin  may  occur 
in  localized  regions  of  accretion  disks  that  perturb  the  radiation 
characteristics  from  the  standard  model  average.  lUE  spectra 
of  UX  Ursae  Majoris  show  substantial  time-dependent  SED 
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changes  for  that  NL  system  (Linnell  et  al.  2008b).  The  adherence 
of  V3885  Sgr  to  the  standard  model  in  this  study  may  indicate  an 
infrequently  violated  average.  On  the  other  hand,  the  observation 
times  of  nonconforming  systems  may  occur  at  times  when 
transients  are  prevalent  but  detailed  frequent  monitoring  might 
discover  times  when  their  SEDs  fit  the  standard  model. 

7.  SUMMARY 

We  use  the  Knigge  (2006,  2007)  period-secondary  mass 
relation  to  assign  the  secondary  mass  and  the  mean  WD 
mass  from  the  same  paper,  consistent  with  an  observationally 
determined  mass  ratio,  to  set  the  V3885  Sgr  WD  mass.  Based 
on  a  analysis  using  38  system  models  and  corresponding 
synthetic  spectra,  we  find  that  a  standard  model  accretion  disk 
with  M  =  5.0  zb  2.0  x  IO-^Mq  yr-\  a  57,000  zb  5000  K 
0.7  Mq  WD,  and  an  orbital  inclination  of  65°  zb  2°  provide  a 
system  synthetic  spectrum  which  optimally  fits  observed  FUSE 
and  STIS  spectra  of  V3885  Sgr.  The  result  is  robust  for  small 
changes  in  the  assumed  M^d-  A  Hipparcos  parallax  fixes  the 
divisor  enabling  the  system  synthetic  spectrum  to  be  matched  to 
the  observed  spectra  and  establishes  the  fit  on  an  absolute  fiux 
basis.  The  model  thus  provides  realistic  constraints  on  M  and 
Teff  for  a  large  M  system  above  the  period  gap. 
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